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A distorted-wave version of the renormalisation group is applied to scattering by an inverse- 
square potential and to three-body systems. In attractive three-body systems, the short-distance 
wave function satisfies a Schrodinger equation with an attractive inverse-square potential, as shown 
by Efimov. The resulting oscillatory behaviour controls the renormalisation of the three-body inter- 
actions, with the renormalisation-group flow tending to a limit cycle as the cut-ofi^ is lowered. The 
.^jl ' approach used here leads to single- valued potentials with discontinuities as the bound states are cut 

f^ . off. The perturbations around the cycle start with a marginal term whose effect is simply to change 

t"*^ ' the phase of the short-distance oscillations, or the self-adjoint extension of the singular Hamiltonian. 

rvj ' The full power counting in terms of the energy and two-body scattering length is constructed for 

_i , short-range three-body forces. 
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I. INTRODUCTION 

The successful application of effective field theories (EFT's) to two-bodyscattering has revived interest in developing 
I> ' model-independent treatments of few-body systems, as reviewed in Refs. [llQl- In the case of two-body systems EFT's 
OO . provide a framework within which the old idea of the effective-range expansion [3|, |j| can be extended systematically 
to describe electromagnetic or weak couplings. Their application to three-body systems requires the addition of 
three-body forces. The observation that such forces are essential for describing low-energy three-body observables 
goes back to the work of Phillips [5J. One model- independent way of introducing them is the boundary-condition 
^^ \ method developed by Brayshaw Q. The EFT approach provides a more practical way of doing so, and one which can 
f^ ' be extended to include couplings to external currents. 

"^ [ In the case of two-body scattering, these EFT's have been well explored. If there is a clear separation between the 

-4— > . low-energy scales of interest and the high-energy scales characterising the underlying physics, then the interactions 
,_^ ' terms in the effective Lagrangian can be organised systematically as an expansion in powers of ratios of low-energy 
^ , to high-energy scales. For weakly interacting systems the power counting in this expansion is just that of naive 
r^ ' dimensional analysis. Since this is the same as the counting in chiral perturbation theory [7J,|3, it is usually known 
.. \ as "Weinberg counting" . 

^ ■ In contrast, in strongly interacting systems with shallow resonances or bound states, simple dimensional analysis 

k>( \ is no longer appropriate because of the appearance of new low-ener gy s cales. These result in a need to resum certain 

'V^ . terms [3, 13 • This has been done within various frameworks ^ ^ly, [13, [l3 1 ^ill of which lead to the the same power 

Cd ' counting. The resulting scheme is often referred to as "KSW counting" . For systems with short-range forces only, it is 

in fact equivalent to the effective-range expansion developed by Bethe and others |3,I3- This expansion is the relevant 

one for few-nucleon systems at low energies, where the deuteron and the ^5'o virtual state lie very close to threshold. 

It may also apply to atomic systems where Feshbach resonances can be tuned to give large scattering lengths. 

Both of these power-counting schemes can be understood using the Wilsonian renormalisation group (RG) [ll|. 
Fixed points of the RG describe scale-free systems. The RG flow near these points can then be used to define a 
power counting for organising the perturbations around them. Weinberg counting arises from the expansion around 
the trivial fixed point of the RG, corresponding to a noninteracting system. KSW counting, or the effective-range 
expansion, is associated with a second, nontrivial fixed point corresponding to a system with a bound-state exactly 
at threshold. 

The RG approach can also be extended to systems with known long-range interactions, such as the Coulomb or 
one-pion-exchange potentials, provided one has identified all the low-energy scales. In Ref. [lj| we applied this to 
examples including Coulomb and repulsive inverse-square potentials, for which well-defined distorted waves (DW's) 
exist. We refer to this extended version as the distorted- wave RG (DWRG). In it, a cut-off is applied to the basis 
of DW's of the known long-range potential so that that potential is unaffected by the cut-off. The advantage of the 
method is that it provides a clean separation between the short- and long-range physics. The resulting nontrivial 
fixed point corresponds to a DW or "modified" version of the effective-range expansion |3i[l3- I'^ this paper, we use 
the DWRG to determine the power-counting for three-body forces. A very brief account of these ideas was previously 



presented in Ref. [ig and a much more extensive one can be found in the Ph.D. thesis, Ref. [T7| . 

The EFT description of systems of three particles must include three-body forces, represented by six-point interac- 
tions in the Lagrangian. To develop a model-independent treatment of these systems, we use a two-body EFT with 
couplings fixed by two-body data, and augment it with three-body terms. However a consistent power counting needs 
to be determined for these new interactions. 

Complications arise in three-body systems because there can be long-range forces resulting from the exchange of one 
of the three particles. The range of these forces is controlled by the two-body scattering lengths. In weakly interacting 
systems, where the size of the scattering lengths is determined by the scale of the underlying physics, these forces can 
be treated as short-range and Weinberg counting applies to the three-body interactions. On the other hand, systems 
with shallow bound or virtual states have large two-body scattering lengths and hence long-range forces. Hence for 
systems where the two-body forces can be organised with KSW counting, the correspnding power-counting for the 
three-body forces is difficult to determine. 

Bedaque et al. [13, [l3 were the first to look at the problem of short-range three-body forces in EFT's. Their 
method is based on the equation originally derived by Skorniakov and Ter-Martirosian (STM) |23| for systems with 
zero-range two-body forces. In higher partial waves and in s-wave nuclear systems with spin or isospin | they found 
that the three-body interactions are irrelevant (in the technical RG sense that they vanish as the cut-off is lowered) . 
This is because centrifugal forces or the Pauli exclusion principle act to ensure that these systems are insensitive to 
short-range physics. In contrast, for three bosons in s-waves or three s-wave nucleons with total angular momentum 
and isospin ^, there is nothing to prevent all three particles from coming together. Bedaque et al. found that three- 
bod y fo rces are not merely important in these systems, they are in fact essential for producing well-defined results 
|18lll9l | (see also Refs. [23,|21|). More recently, a power-counting scheme for the three-body forces in these systems 
has been constructed by renormalising the STM equation order-by-order in the energy |23 | . Our work confirms this 
power counting within the framework of a full RG analysis. 

In other recent work on this problem, Phillips and Afnan [2Jl have looked at the STM equation in the case of three 
s-wave bosons. They were able to reproduce the results of Bedaque et al. by using a subtractive renormalisation and 
introducing a single piece of three-body data, namely the three-body scattering length. The introduction of this piece 
of data is equivalent to the LO three-body force of Ref. jl£| . Similar results have also been obtained by Mohr '^ . 

An alternative to the STM equation for three-body systems with contact interations is provided by the work of 
Efimov ^d| (for a review, see Ref. [231). In general, Efimov's equations are difficult to solve (see, for example, |23]) 
because they involve a nonseparable boundary condition. However, in the case of infinite two-body scattering length, 
the equations become separable in suitable cooordinates and can be solved exactly. In this limit the Hamiltonian 
becomes scale free, which should not be surprising given that this is the nontrivial fixed point of the two-body system. 
In hyperspherical coordinates, the corresponding three-body Hamiltonian contains an inverse-square potential (ISP). 

The strength of the ISP is determined by the statistics of the system. In systems with nonzero orbital angular 
momentum or where the Pauli exclusion principle keeps the particles apart, the potential is repulsive. However if 
there is no angular momentum or exclusion principle, as for three s-wave bosons, there can be an attractive ISP. This 
corresponds to the cases where Bedaque et al. found a need for a three-body force. 

Because of its relevance to three-body EFTs, the attractive ISP has been the subject of several recent papers |29ll30l 
l3l| . It will also play a central role in this paper, since the DWRG analysis for scattering in the presence of this potential 
will also provide the basic power counting for the three-body system. The attractive ISP is particularly interesting from 
a mathematical point of view because it is sufficiently singular that its wave functions have a logarithmic oscillatory 
behaviour near the origin |32| . This makes it impossible to define a "regular" boundary condition on the wavefunctions 
at the origin. Mathematica lly, the resolution of this quandary is well known: we need to form a self-adjoint extension 
of the Hamiltonian |33l l34l l35l | . More physically, one can think of this as introducing a boundary condition to fix the 
phase of the oscillatory solutions near the origin. 

Although choosing a self-adjoint extension ensures that all observables are uniquely defined, it does not lead to a 
ground state. Instead the bound states form an infinite tower with geometrically spaced energies. This pattern was 
found by Efimov in the context of three-body systems po], although the lack of a ground state was noted first by 
Thomas 36J. The pattern is not quite scale-free since it requires the input of a new scale, such as the energy of one 
of the bound states. This scale is provided by the self-adjoint extension. 

An alternative to introducing a boundary condition to fix the self-adjoint extension of the ISP is to regulate the 
singularity of the potential at short distances and introduce a counterterm 29, 37J. In essence this is what the 
momentum-space regularisation of the STM equation does for the three-body problem [l9ll2lll22l. A gain, one piece 
of three-body data is needed to fix the scale in the extension, as also noted by Phillips and Afnan [2J| . Other authors 
have shown explicitly how a short-range force can be used to select an extension 29, jsj, |33 , but have not developed 
the power counting for all possible terms in that force. Furthermore, the three-body forces obtained in this manner 
turn out to be multi- valued, with no obvious procedure to resolve this |29l l3ll| . 

An interesting feature of the approaches based on regularisation of the singular interaction is that they exhibit a 



cyclic behaviour as the scale of the cut-off is varied p9, 21, 2</,!33- In fact, as pointed out by Wilson ^8^, they form 
a novel kind of limit cycle of the RG. Several examples of this have only recently been highlighted ,39. .4(1 l4ll| . In 
addition, Braaten and Hammer have pointed out that a rather minor fine tuning of the quark masses in QCD would 
bring the infrared limit of the theory to the effective-range fixed point in the two-nucleon channel, and hence to a 
limit cycle for three nucleons |43| . 

In this paper we apply the DWRG method, first to the attractive ISP corresponding to three-body systems with 
infinite two-body scattering length, and then to more general attractive three-body systems. We apply a cut-off to the 
deeply bound states as well as the continuum, and this leads to an RG equation with single-valued but discontinuous 
solutions. As already mentioned, the results confirm the power counting found by Bedaque et al. |23|. By taking 
advantage of Efimov's separation of the full three-body wave function, we are able to show that this power counting 
is not limited to the STM "slice" through it, where two of the particles coincide. Because of the clean separation of 
short- and long-range physics in the DWRG, we are able to derive this counting in a more transparent way. 

The lowest-order three-body interaction is marginal, in the sense that its RG flow is only logarithmic in the cut- 
off. As we shall explicitly show, it represents the same degree of freedom as the choice of self-adjoint extension. The 
resulting three-body interaction is built out of perturbations around a single- valued limit-cycle solution of the DWRG. 
By constructing the corresponding physical scattering amplitudes, we are able to directly relate the coefficients in the 
three-body potential to observables, through a version of the DW effective-range expansion. 

II. THE RG FOR THREE-BODY FORCES 

For simplicity, we consider here a system of three bosons of mass M , in a state of zero orbital angular momentum. 
We assume that the two-body interaction produces a single shallow bound state with binding energy 7^/M , corre- 
sponding to a "binding momentum" i^. We start by considering the effects of two-body forces only and we denote the 
corresponding three-body Green's function by G2 (j>) , where p is defined in terms of the total centre-of-mass energy 
hy E — p^ /M. This function has the spectral decomposition. 
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where ^(x) is the unit step function. Note that in a system with two-body contact interactions, we require an 
additional boundary condition to define the DWs and hence the Green's function, as discussed below. 

Three types of wave appear in this decomposition. In the first term, we have the bound states of three particles, 
|\E'„). In the final term, we have states with three incoming free particles, l^q.fc), which we label by their total centre- 
of-mass energy, E = q^ /M , and relative momentum, k. In the middle are states with one incoming free particle and 
a bound pair \^q^i^). These states are normalised as described in Appendix A. 

In this paper we work with the coordinate-space representation of the DW's. Our EFT's are expressed in terms 
of contact interactions and so we are interested in the region where all three particles are close together. A natural 
measure of the proximity of the three particles is the hyperradius, R, defined in Appendix A. We write the DW's 
as functions of R and the five hyperspherical angles, collectively denoted by i7, which represent the other degrees of 
freedom in the centre-of-mass frame. The precise specification of these angles is not needed since all f2 dependence 
will factor out of our results. 

Now consider the effect of introducing a three-body interaction, V3. This could be inserted directly into the Faddeev 
equations [43. However, for our purposes, it is sufficient to use the "two-potential trick" |22| to define a T- matrix for 
the additional scattering produced by V3. We write a full 3 — » 3 T- matrix, T(p), in terms of the one with two-body 
forces only, T2{p), plus an additional piece, T^{p), that acts between the DW's for the two-body forces: 



{q,k\T{p)\q\k') = ((?,fc|r2(p)|g',fc') + (*,,fe|f3(p)|vl',,,fc, 
The DW term T^{p) satisfies the Lippmann-Schwinger (LS) equation, 

f:i{p) = Vz + V^G2{p)f:i{p). 



(2) 



(3) 



There is no problem with connectedness in this equation since the kernel contains the three-body force. 

For zero-range two-body interactions, the small- i? behaviour of the wave functions in Eq. (^ can be found from 
the approach of Efimov |26| . The essential elements of this, in our notation, are outlined in Appendix A. Since the 
wave functions do not have well-defined limits as i? ^ 0, we choose our effective three-body interaction to act at some 
small, but nonzero hyperradius, R ~ R. At this point the separations of all three particles are less than or of the 



order of R. This additional regularisation is quite separate from the running cut-off which wiU give the RG flow. As 
with the singular long-range potentials studied in Ref. |l4|. it is needed to make all wave functions appearing in the 
RG equation well-defined. Provided we choose R to be in the region where the waves have reached their asymptotic 
small- i? forms, the common dependence on R can be factored out of the problem. One could therefore imagine taking 
the i? — > limit of the results, corresponding to a zero range force. In addition, we take the interaction to act at 
some fixed set of hyperangles il = il. Since the wave functions at small R are separable, all hyperangular behaviour 
can be factored out and this arbitrary choice will not affect the results. We can therefore define 

{^pMV3\-^p,k') ^ R^%,kiR,^)'^P,k'iR,n)V3{p,k,k'). (4) 

In an EFT we aim to integrate out all the unknown short-range physics and replace it by a three-body force which 
can be expanded in powers of the low-energy scales of the system. The tool which allows us to determine the scaling 
behaviour of the terms in V3 is the renormalisation group (RG) [llj . The first step in setting this up is to impose a 
cut-off. A, to separate the low-energy states which we wish to treat explicitly from the high-energy states which will 
be integrated out. We then renormalise the theory by making the effective interaction A-dependent and demanding 
that observables should be independent of the cut-off. 

As in the approach developed in Ref. [Ij|, we impose this cut-off on the DW's of our pairwise forces. By applying 
the cut-off to these waves rather than the apparently simpler free waves, we ensure that the role of the three-body force 
is related purely to short-range three-body physics. If the cut-off were applied to the free waves, it would also remove 
parts of the physics associated with the two-body forces. This would have to be compensated by the three-body force, 
which would then satisfy a far more complicated evolution equation. 

We impose the renormalisation condition that amplitude T^ be independent of cut-off, 9a Ts — 0. Differentiating 
the LS equation Q for T^ then leads, after the elimination of Ts, to a differential equation for V^, 

The RG equation is then obtained by rescaling the potential V3 and and all the low-energy scales in this equation. As 
described below, the definition of the rescaled potential V3 depends upon the form of the DW's close to the origin. The 
boundary conditions on V3 are that it should been analytic in all rescaled energies and momenta. These follow from 
our requirement that the terms in the three-body force arise from local six-point vertex terms in the EFT Lagrangian. 
The A-dependence of the corresponding terms in V3 then allows us to classify them according to some power counting. 

So far we have said nothing about the nature of the two-body interaction, except that it leads to a single low- 
energy bound state. This could be represented as a contact interaction, whose terms correspond to the effective-range 
expansion. Equivalently it can be thought of as a boundary condition on the logarithmic derivative of the wave 
function when two of the particles coincide. Combining this with the overall symmetries of the wave function leads 
to Efimov's boundary condition 26], as outlined in Appendix A. 

This boundary condition becomes separable in the limit where the hyperradius, R, is much smaller than the two- 
body scattering length, 02. The Hamiltonian is then diagonal in the hyperangular momentum variable, s. The 
potential in each hyperangular- momentum channel is an ISP (as it must, since all scales have been eliminated from 
the problem) . The resulting Schrodinger equation can be solved analytically in terms of Bessel functions of general 
order. 

In the cases of three bosons in s- waves or two neutrons and a proton with total spin ^ , the lowest hyperangular 
momentum gives rise to an attractive ISP. As discussed in the introduction, there is an ambiguity in the wave functions 
for this potential. In order to resolve this, we need to introduce an additional boundary condition to form a self- adjoint 
extension. As shown in appendix A, the form of the resulting DW's at small R is 

^{R, n) ^ 60 i^) sin(so ln(p,i?) - 0) , (6) 

where Sg ~ 1.006 is the magnitude of the smallest hyperangular momentum, 6 is defined in Ea. ljA14p . ^S(j(r2) is the 
hyperangular wavefunction, and p* is the scale introduced by forming a self-adjoint extension. This form is invariant 
(up to an overall sign) under the replacement, 

p*->p*e"'^/"°, neZ. (7) 

As a consequence all observables are invariant under this transformation and so all of the values of p* in Eq. Q define 
the same self-adjoint extension. 



III. INFINITE SCATTERING LENGTH 

In the limit of infinite two-body scattering length, the Haniiltonian due to pairwise forces is separable. If we assume 
also that the three-body force is diagonal in the hyperangular momentum then the RG equations for the various 
channels are uncoupled. This approximation neglects parts of the interaction that couple the channels. We introduce 
it here so that we can focus on the lowest hyperangular-momentum channel, which contains the leading interactions. 
Later, when we derive general results, we shall not require it. 

The higher hyperangular-momentum channels, Si with i > 0, correspond to repulsive ISP's. Their RG analysis is 
thus identical to that described in Ref. [lj|- It leads to a power-counting scheme in which the LO term in the force 
scales with (p/Ao)^'*' , where Aq is the scale of the underlying physics. The most important of the interactions based 
upon this power counting is the LO term in the si = 4 channel, but even this term is heavily suppressed. Similarly, 
in systems which do not lead to attractive ISP's, the lowest-order three-body forces are always suppressed by powers 
of p/Aq. Here we concentrate on the So channel in attractive systems, where the effects of three-body forces are 
enhanced. 

A. RG equation 

The hyperradial Green's function in the So channel is given by 

g2{p;R,R) = — >^ -^-3 + -j / dq -, — -^— — . (8) 
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We use lower case to denote quantities in this channel, for example this Green's function and the potential W3, to 
differentiate them from those for the full problem. The DW's Us^ (p, R) are simply solutions of the ISP Schrodinger 
equation l|A9|l . 

The bound states Us^ {R) have energies which, from Eq. (jA17|) . form a tower with geometric spacing, 

E^ = -^e'^^''\ (9) 

These respect the symmetry of Eq. ((JJ. They accumulate at zero energy and extend deeper and deeper with no 
ground state. The shallow states are referred to as Efimov states, although the lack of a ground state had been noted 
much earlier. In 1935, Thomas 3&| pointed this out for three-body systems bound by contact two-body forces. 

Before we can write down the DWRG equation for the three-body force in this channel we must first decide how 
to handle the bound states. In the previous examples studied with the DWRG method [lj|, the cut-off was applied 
purely to the high-energy continuum states. In those examples, any bound states were shallow, with typical momenta 
much smaller than the scale of the underlying physics, and so lay within the domain of the EFT. Here, however, 
the lack of a ground state means that this is no longer true. The deeply bound states are consequences of our use 
of contact interactions to represent the two-body forces, and hence are unphysical artefacts. In keeping in the EFT 
philosophy these states should be truncated and their effects absorbed into the effective three-body force. We chose 
to cut-off the Green's function by removing all states with energy outside the range —A?/M < E < A^/M: 

f R R'^ ^^ V- u^::\R)ui':\R') , M f^^ u,,{q,R)u,,{q,R') ^n^ 

\p\<A P + ^" ^ "'0 p^ - q^ + le 

We shall see that truncation of bound states will allow us to get an RG equation with single- valued solutions. 

Using the separable form, Eq. 0, of the wave functions, we can project the differential equation for the potential 
© onto the Sq channel. We can then absorb the common angle-dependent factor |^so(^)P i^ito into the potential by 
defining 

vsip,A)^\^,MfV3ip,A). (11) 

Note that we have chosen to focus on potentials that depend only on energy since, as discussed in Refs. [l2 . ll4l |. these 
contain the leading perturbations. 

Projecting the differential equation jSJ onto this channel and inserting Eq. (|10|l we obtain an equation for V3, 
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«3b,A)2. (12) 



The first term is produced by cutting off tiie continuum states, and is similar to the expressions in Refs. [l^iliJI- The 
series of discontinuities represented by J-functions results from the truncation of the bound states. 

Finally, to obtain an RG equation we need to rescale any low-energy scales in the problem. In this case we have 
only one, the on-shell momentum p, and we define p = p/A. To form a dimensionless potential, we multiply by 
the mass M. We also take the expressions for the wave functions near the origin, Eqs. I|A13I fSTOI) . and absorb the 
i?-dependence into the potential by defining 

?)3(p,A) ^ sin2(soln(p*i?)-0)w3(pA,A). (13) 



The rescaled potential v^ satisfies the RG equation 
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where p„(A) ~ pn/A — p*e'"^/'*''/A. Note that we write p„(A) as a function of A since we do not rescale the scale 
from the extenstion, p^, and so pn varies with A. As noted before |l4l |. this type of equation is more conveniently 
rewritten as a linear equation for l/Ds, 

9/l\ ^ d / 1\ sinh(7rSo) -k ■^ 1 xt- f\\ \ 

^dl [~J ^^dp [¥,) ~ [cosh(^So)-cos(2s„ln(A/p,))l(l-p2) +;; 2- TT^'^^P"(^^-l^ (^^^ 
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This RG equation is the same as that for an attractive ISP in two dimensions. In fact an attractive ISP in any 
number of dimensions leads to an equation of this form, because the different real power of R which appears in the 
short-distance wave functions exactly compensates for the radial factor in the Jacobian. All that differs is a numerical 
factor in the definition of V3, associated with the angular integration. 

B. Solutions 

The physically acceptable short-range potentials are given by solutions to the RG equation (|14|l which satisfy the 
boundary condition of analyticity in p^ for small p. We look first for fixed points, A-independent solutions of the 
equation. The power counting for terms in the potential can be determined from the perturbations around a fixed 
point that scale with definite powers of A |13. [lj|- These are eigenfunctions of the linearised version of the RG 
equation. 

An obvious fixed point is the trivial one, v^ ~ 0. The power counting based on it can be found by substituting 

V3 = CA^0(p) (16) 

into Eq. H14|l . and linearising to get the eigenvalue equation 

P^-.0. (17) 

dp 

This equation is easily solved, giving (/){?) = P^- Imposing the boundary condition of analyticity in jp, leads to the 
eigenvalues /x = 0,2,4,.... Hence the general solution in the region of this fixed point is 

00 

*3(A,p) = ^C2„A2V"- (18) 

n=0 

The leading term in this expansion is marginal, that is, it does not scale with any power of A. We therefore expect to 
find logarithmic dependence on A associated with this perturbation. To resum these logarithms we need to construct 
solutions of the full nonlinear RG equation. 

In fact the presence of A-dependence on the right-hand side of Eq. H14|) means that no other fixed-point solution can 
be found. However that dependence is only logarithmic and so we may look for slowly evolving solutions which also 
depend logarithmically on A. Perturbations about such a solution can still be used to construct a power counting. 

As in Ref. |l4| , our starting point for constructing a nontrivial solution to the RG is the basic loop integral, which 
in this case is given by 

Iip,A)=V dq '^ I "> (19) 

Jo |_cosh(7rso) - cos(2so ln(Aq/p^))\ (p^ - q-^) 
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FIG. 1: The contour C in the complex g-plane used to construct the solution J(p, A) of the DWRG equation. The bound-state 
poles of the integrand occur at g = ip„ = ip»e"'^' ^''/A and the propagator poles at g = ±p. 



By substituting this for l/v^ in Eq. p5|l . it is straightforward to verify that it satisfies the continuous version of the 
equation (without the discontinuities of the final term). However it is not yet an acceptable solution since, apart from 
missing the bound state discontinuities, it is nonanalyitic for small p. 

The propagator pole at q = p in the integrand approaches the endpoint of the integral as p — > 0. This leads to loga- 
rithmic dependence of / on p. To avoid this troublesome endpoint, and also to introduce the required discontinuities, 
we need to find a similar integral over some contour in the complex g-plane avoiding the singular region. 

The first step is to rewrite the integrand of / as 
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This integral has the same nonanalytic behaviour on the real axis as that in Eq. H19() . However, since it is no longer 
associated with an endpoint of the integral, we are free to deform the contour of integration into the complex plane 
to avoid the dangerous region around q = 0. We therefore define 



j(p,A)= / dqhiq), (23) 

Jc 

where C is the contour of integration shown in Fig.^ This runs from g = — 1 to ^ = 1 in the upper half plane and 
crosses the imaginary axis aA, q = i. 

Fig. n also shows the pole structure of the integrand h{q). Apart from the two propagator poles at q — ztp, there 
are bound-state poles at 



<? = zp„(A)=*^e"^/^°. 



(24) 



These poles are important since their positions vary with A and so they can cross our integration contour. When 
this happens they generate discontinuities in j(p. A) . By choosing our contour to cross the imaginary axis at q = i 
(corresponding to an energy E = — A^/M) we can arrange for these discontinuities to match precisely the ones we 



require in Eq. (|15|l . To see this, note that as A is varied from p„ — e to p„ + e, the pole at g = zp„ crosses the contour 
and produces a discontinuity 



A=p„+e 7r 

= -2Trin[h{q),ipn] ^ — 

A=p„-e So 1 +p"^ 



i(p,A) , =_27ri7^[/l(<z),^p„]--^^2, (25) 



where TZ[f{z),ZQ] denotes the residue of f{z) at the pole zq. This is precisely the strength of the (5-function at 
Pn{A) = 1 in Eq. ^. 

The integral j{p, A) satisfies the full RG equation (|15|) and is analytic about jP = 0. We therefore define the rescaled 
potential by 

4°\p,A)^[j{p,A)]~\ (26) 

This is not a fixed point since it does depend logarithmically on A, but in it these logarithms have been resummcd 
to all orders. The invariance of the system under Eq. Q means that this dependence on In A is periodic. As pointed 
out by Wilson 38], this periodic behaviour provides an example of a limit cycle of the RG. 

Like a fixed point, the limit-cycle solution can be used to define a power counting for the perturbations around it. 
Because the RG equation in the form of Eq. H15|) is linear in l/ua, it is possible to construct exact solutions containing 
these perturbations. 



hip, A) 4°^(p,A) 



^ oo 

+ ^C2„A2V". (27) 



r!=0 



The perturbations around the cycle thus have exactly the same power counting as those around the trivial fixed point. 
The leading, energy- independent term is marginal. Each additional power of energy (p^) gives a term two orders 
higher in A. 

The general solution, Eq. (|27|l , shows that there is a family of limit cycles parameterised by the marginal perturbation 
Cq. This lack of uniqueness arises because the marginal perturbation is A-independent and satisfies the homogeneous 
version of the RG equation H15|l . As a result, an arbitrary amount of it can be added to the limit-cycle solution, 
Eq. H26() . The leading perturbation around any cycle is marginal because it pushes the solution into a nearby cycle. 
In the limit Co — + oo, the entire cycle is compressed into the trivial fixed point v^ = 0, which can be regarded as the 
limiting member of the family of cycles. All other perturbations vanish as A — > and hence are stable. Any more 
general solution will thus tend to one of the limit cycles as A — > 0. 

Each limit cycle contains one discontinuity every period, when a bound state is removed from the low-energy 
domain. The point in the cycle where this occurs depends on how we truncate the sum over bound states. Here 
we have chosen a cut-off that is symmetric between positive and negative energies, although other prescriptions are 
equally acceptable since physical observables should not depend on them. 

The truncation of bound states ensures that our solution of the RG equation is single- valued. Our choice corresponds 
to a particular path from g = —1 to +1 for the contour C used to evaluate j{p,A). Other, topologically different 
paths would correspond to different branches of a general multi-valued solution. Such multi-valuedness appears in 
other methods used to renormahse the three-body potential, as in Refs. [13. 12(1 I2lll2 ^ . 

C. Scattering observables 

To relate the parameters appearing in our solution, Eq. H27|l . to scattering observables, we need to construct the 
corresponding T-matrix. If we use the separable form of the wave functions, we can define a DW T-matrix t^ in an 
analogous manner to Eq. I|12|) . Projecting the LS equation Eq. (PJ onto the sq channel, we get the correspondimng 
equation for ^3. With our choice of energy-dependent potential, this can be solved directly to get the on-shell T-matrix, 

{u;\k{p)\up ^ e'''^-^P^\up{R)\'v3{p,A) [l~g2{p,A;R,R)v3{pA)r\ (28) 

where the superscripts ± denote waves with incoming or outgoing boundary conditions, and 62 is the phase shift in 
the So channel produced by the pairwise forces alone. 

Unitarity and conservation of hyperangular momentum allow us to express this in terms of a phase-shift as 

„2i<52(p) />/- 

= -— (cot<53(p)-*). (29) 



(up \t3{p)\uj) 2ttp 



Here ^3 is the additional phase shift produced by the three-body force. It is related to the full phase shift 5 by 
$3 =5 — 62- Using the fact that our potential acts aX R = R, we can combine Eas. (|28l I29|) to obtain an equation 
relating V3 and the phase shift: 

A'l 
\up{R)\^ — {cot ~53~i)^g2{p, f^;R,R)-v3{p, A)-\ (30) 

When we use the expression in Eq. (|l()(l for the regularised Green's function and the explicit forms of the wave 
functions from Appendix A, we find that the R dependence can be factored out to leave 

n I sinh(7rSo) \ / ? , s \ f gsinh(7rsQ)(i(7 

{cotS3{p)-t) - 



2 \ cosh(7rso) ~ cos(2s(,ln(p/p*)) / 7o [cosh{T:so) ~ cos[2so\n{q/p^))]{p'^ — q'^ + ie) 

+- E -^-^-Anrv (31) 

^OiT^aP+P" v3ipA,A) 

To write this in a form which is explicitly independent of A, we introduce the rescaled variable q = q/A and substitute 
in our solution for v^ from Eq. (|27() . This gives 



/I !• CO 

dqh{q)- dqh{q)-27Ti ^ 7^ [/l(g), ip„/A] - ^ C2„p'", 
'^ '^ Id„I<A n=0 



n I sinh(7rSo) 

2 \^cosh(7rso) -cos(2soln(p/p*))^ „_, „^ 

(32) 
where h{q) is defined above in Eq. H21|) and we used the fact that the sum over the bound states can be expressed in 
terms of the residues at the bound-state poles (see Eq. (|23). 

The two integrals over h{q) can be combined to form a single integral round a closed contour. The ie prescription 
in the first integral implies that the contour of integration goes above the propagator pole at q = —p and below the 
one at q — p. The contour thus encloses the pole at positive q. It also encloses the bound state poles at q = ipn/A 
with pn < A. Using Cauchy's theorem, we find that the residues from these exactly cancel the sum over bound states 
in Eq. ( P|l . 

We are then left with a DW effective-range expansion, 



cosh(7rso) — cos(2?7(p)) I " '" ' \ ' i: / tt 



00 



' (cot Saip) -i)= cot 7j{p) + ^ - - E ^2nP'", (33) 



n=0 



where r]{p) — —Sq ln(p/p*) is defined in Eq. I|A15(I . Make use of trigonometric additon fomulae, one can show that the 
imaginary parts of the left and right sides are equal. In this expansion, all nonanalytic behaviour has been subtracted 
or factored out into the trigonometric functions of In p. The remaining energy dependence can be expanded in powers 
of p^ and the terms correspond directly to perturbations around the limit cycle. 

The roles of the limit-cycle potential and its marginal perturbation Cq are not obvious from Eq. H33|l and so, to 
clarify these, we construct the total phase shift, 6. Making use of Eqs. H33|l and (|A16|I . the corresponding full ^-matrix 
can be written 

Z{p) ' ^^ ^ 

where 

Zip) = cos Up) + '-^)^l sin Up) + ^) f: C2np'-. (35) 

For the limit-cycle solution of Eq. (|26|) with no perturbations (all 6*2,1 = 0), we have 

Z{p) = cos (^vip) + ^) . (36) 

Comparing this to Eq. (jA16|) for the pure long-range force, we see that the phase ijip) of the waves near the origin 
has been shifted by tt/2. This implies, for example, that the bound states lie at p = ie~'"^^^°Pn, which correspond to 
the geometric means of the bound state energies for U3 = 0. This potential thus corresponds to the system which is 
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"furthest away" from the one with ws = 0, in the sense that it leads to the maximum possible changes to physical 
observables. 

To elucidate the role of the marginal perturbation, it is convenient to express its coefficient in terms of an angle a^ 

Co = -|cotCT, (37) 

and to redefine the other short-distance parameters according to 

2 V^ ^ 2n esc a X)n=l ^2nP " (oo\ 

TT ^^ sm cr + COS a Y.n=l CLp^ 

In terms of the new parameters, a and C2„ , n > 1 we find that the S'-matrix can still be written in the form of 
Eq. ijnil), but with Z{p) replaced by 

Z\p) = sin Up) + a + ^) + cos Lp) + a + ^) £ 0,^^. (39) 



n=l 



This form makes it clear that a (or rather Co) just has the effect of shifting the phase rjip). Scattering observables 
depend only on the combination 

pi = e-'^/^>„ (40) 

and not on p^ and Cq separately. This shows that p, and Cq play the same role in determining the phase of the 
wave functions for i? — > (the self-adjoint extension of the Hamiltonian). We can use Cq to change the self-adjoint 
extension from the initial one specified by p* to any other. In particular, a — (corresponding to the trivial fixed 
point) leaves the initial extension unchanged, whereas a = tt/2 (Co = 0) produces the largest possible change in the 
extension, as already noted. Furthermore, there is a one-to-one mapping between all possible limit-cycle solutions 
(obtained by varying Cq from — cx3 to -|-oo or equivalently a between and tt) and the self-adjoint extensions. This 
is in marked contrast to the relationship between p* and the self-adjoint extensions, where infinitely many equivalent 
p*'s give the same extension. 

The bound states of the system are given by the zeros of Z'{p). For any short-range potential, the bound states 
still accumulate at zero energy, as expected since they are consequences of the inverse-square tail of the long-range 
potential. The shallower bound states are insensitive to the short-range perturbations, C2„ with n > 1. Their positions 
are controlled by p'^ rather than the original p*. If we choose p^ to give the correct shallow Efimov states, we can set 
cr = and use the expansion around the trivial fixed point. This corresponds to expanding a DW iiT-inatrix and has 
the form 

tan^3(p) ^ E^^i CkP^" ,^^. 

sinh(^so) cosh(^So)-cos(277(p))+sin(2r/(p))^^^,C^„p2»' ^ ^ 

The deeper bound states are of course strongly affected by the higher-order perturbations and will not in general 
follow the simple constant-ratio pattern of the Efimov states. At some point they will fall outside the range of validity 
of the EFT and so we can say nothing about them, except that some short-range physics must act to ensure the 
existence of a ground state. 

IV. FINITE SCATTERING LENGTH 

In the general case of finite scattering length it no longer makes sense to expand observables in terms of the 
hyperangular wavefunctions, instead we must consider a much more general RG equation. In particular, the low- 
energy scales now include 7 — 1/02 where 02 is the two body scattering length. This can also be thought of as the 
typical momentum in the low-lying bound or virtual state of two particles. The coupling between waves with different 
hyperangular momenta means that the states should now depend on the relative momentum fc of a pair, as well as 
the total energy, expressed as a momentum p. Nonetheless Efimov's separation of variables still applies whenever the 
distances between the three particles are all much less than the two-body scattering length [2g. It therefore controls 
the short-distance behaviour of the three-body wave function in an EFT with contact interactions and hence the 
scaling of the short-range three-body interactions. For attractive three-body systems, the the basic structure of the 
RG evolution will be similar to that in the previous section since the wave functions at small R will still be dominated 
by the Sq channel. 
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The truncated Green's function is given by Eq. |^ restricted to states inside the energy range —h?/M < E < A? /M: 



M ^ |vI/„)(*. 



^^p^A P'+P' 






|'J'g,.7)(*9,nl+^('?')- / dk\^q,k){^q.k 



Inserting this into Eq. (0) and using Eq. Q we obtain a differential equation for the potential 



(42) 



dA 



Mm 



27r2 p2_A2 



V^{p,-|,k,ir,A)\^A,^^{RM''V^{pn,il,k'■,K) 



dk"V3{p,j,k,k";A)\^A,k"{R,m^V3{p,^,k",k';A) 







MR'^ 



E 



^^liI^l^SiA-p^)V3ip,j,k,tpj3;A)V3ip,j,iPn/S,k';A). 



47r ^—^ p2 ^ p 
For i? <C a2 all the DW's will be dominated by the s,, channel and will therefore tend to the forms 

\%ARM' - -VsiP,kn,P.)\m r'''^'"'t*''^-'\ 

p R^ 



i^niRM' - piv^^\i,p.)\m)\ 



i?4 

i2sin^ (soln(p»i?) 
i?4 



-,(") 



(43) 



(44) 



The various normalisation functions, I?3, I?2 and "D^ , are determined by the external boundary conditions on the 
DW's. They can be found by solving the full Faddeev equations (or equivalents). These functions are essential to the 
RG discussion, since it is through them that information about the long-distance physics is communicated to short 
distances. 

The forms of the DW's above were chosen to ensure that the normalisation functions are dimensionless. They have 
a common dependence on the hyperangles given by the function 



m) = Ao E 



2 . , fSolT 

-^ sm(2(/3i) V 2 



(45) 



where ipi is the hyperanglc defined in appendix A. Since this dependence can be factored out of the Green's function 
6*2 (p) at small i?, it is unimportant to the RG discussion. In this context, it is convenient to define a Green's function 
with this common short-distance behaviour divided out. 



^2(P,7,P*) 



27r2i?4 



G2{p;R,n;R,n) 



M \^{Cl)\^sm^{soHp,R)^9)' 
Using the forms of the DW's given in Eq. H44|) we define our rescaled potential by 



(46) 



M 



V3{p,j,k,k';A) = -^lam^sin^SoHp^R) ^ e)V3{p,j,k,k';A), 



2n 



(47) 



where, as usual, p — p/A, etc. We also require rescaled versions of the normalisation functions V, which we define by 

p3(^,7,A)=p3(A,Afc,A7,p,), p2(7,A)=I?2(A,A7,p,), ^^'^(7, A) = I?^')(A7,p,), (48) 

Since p, only occurs in these dimensionless functions as a ratio with A (the only other scale) we have suppressed the 
dependence on it to simplify notation. We also define a similarly rescaled Green's function, 



^2(7, A) =e^2(A,A7,p, 



(49) 
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The rescaled potential V3 satisfies the RG equation 

A^TT- = p^TTT + 7^7:^ + fc — - + k ^^ 
dA dp 97 dk dk' 



1 



1-p^ 



t>3(p, 7, fc, ^7; A)P2(7, A)T>3(p, 7, i7, fc'; A) 



+ - / *"]/3(p,7,^,fc";A)p3(fc",7,A)]/3(p,7,^",fc';A) 



TT 







00 _. 

The boundary conditions on V3 are that it be analytic in p^, 7, k^ and k'^. Again these follow from our requirement 
that the potential arises from six-point contact terms in the effective Lagrangian. 

Note that in the P's and the rescaled bound-state momenta j5„ (A) the only explicit dependence on A occurs in the 
ratio A/p^. (The scales p* and A are the only ones in these dimensionless functions.) The scale p, is introduced, as 
before, by the self-adjoint extension which fixes the phase of the waves for R —^ 0. All quantities appearing in the 
RG equation are thus invariant under the transformation of p* in Eq. and hence are periodic in p* . Furthermore, 
since p* and A always appear as a ratio in the RG equation, that equation must be unchanged under 

A^Ae-"^/'°, neZ. (51) 

In interpreting this for the bound state momenta p„(A) we need to be careful. Depending on how we label the states, 
these momenta may not obey the symmetry individually. However the complete set must always obey it collectively. 

To see how this happens, consider the entire rescaled spectrum, p„, at some initial scale A = Ai. As A decreases the 
rescaled bound state momenta increase, until at A = Aie~'^^^° the spectrum has shifted so that Pn{A) = p„+i(A'). 
This is similar to what happened in the case of infinite scattering length, except that a new bound state must have 
appeared at threshold and moved down to po(Ai). This is possible because in following the RG flow we keep the 
rescaled quantity 7 fixed, not the physical value 7. The spectrum at A == Aie^^/'^° is thus identical to that at Ai. If 
we label the bound states so that po(A) always refers to the shallowest, then each of the Pn{A) will be invariant under 
Eq. ^. 

The general RG equation is complicated because it includes couplings between the different three-body channels. 
As a starting point for solving this equation, we consider a solution which is independent of the two relative momenta, 
k and k'. This simplifies the equation considerably since it allows us to divide through by V3(p, 7; A)^ to obtain a 
linear equation for I/V3, 

4 (^) -4 i^X (4) - T^^<-^' + ||<H.A)T^^(-(A) - ')> (^^) 

where we have introduced the function 

P(7, A) = Ail, A) + - f dk"f>,{k'\ 7, A). (53) 



This function can be thought of as the (rescaled) short-distance part of the projection operator onto continuum states 
with energy E = A^/M. 

Our strategy for solving Eq. (|52|l will mirror that used in the simpler case of infinite scattering length. As before, 
we can immediately write down a solution to the continuous equation using an equivalent to the basic loop integral, 
Eq. H19|) . Again this will not satisfy the analyticity boundary conditions because of the behaviour of the integrand 
at the endpoint g = 0. Furthermore this solution will not generate the discontinuities at the bound states. We must 
therefore find a generalisation of Eq. (|20(l and write the solution as a contour integral that avoids the point q = 0. 

Guided by the results for infinite scattering length, we define an integral along the same contour C as before, 



where 



J(p,7,A)= / dqHiq), (54) 

Jc 



H{q) = -^^^g2{l/q,m- (55) 

IT p^ — q'^ 
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FIG. 2: The contour C used to construct the DWRG solution J(p, 7, A). The bound-state poles of the integrand H{q) occur 
at g = ipn- The branch cut running down the imaginary axis corresponds to the 2 + 1 elastic continuum. This is joined by a 
cut starting at g = corresponding to the three-body continuum. There are also propagator poles a.t q = ±p. 



The contour C is shown in Fig. |21 along with the singularity structure of H. This structure follows straightforwardly 
from the form of G2 and is derived in appendix B. 

In appendix B we also show that the function P(7, A) in Eq. H52() can be rewritten in terms of Q2 as 



P(7,A)--(g2(7,A)-g2(-7,-A) 
n 



(56) 



This result can then be used to show that j{p, 7, A) satisfies the continuous part of Eq. (|52() . In addition, the bound 
state poles in H cross the contour of integration at p„ = 1 and so produce the necessary discontinuities. Finally, the 
contour never approaches the branch points at q = and q = ij and so the integral is analytic in the small scales. 
This allows us to take our nontrivial solution to the DWRG equation to be 



vr{p,l,A)=[J{p,l,A)]- 



(57) 



The periodic dependence on A of all quantities appearing the DWRG equation (|52|l implies that our solution 
Vg^ (p, 7, A) must also be periodic. Like the analogous solution for infinite scattering length, it is therefore a limit- 
cycle of the DWRG. 

As a result, it must have an infinite number of discontinuities as A — )■ 0, which follow a periodic pattern according 
to this symmetry. At first glance this may seem odd since these discontinuities are associated with the truncation of 
bound states and there are only a finite number of bound states in the truncated Green's function. However as A is 
lowered for fixed 7, new bound states appear in the rescaled spectrum as discussed above 

Now that we have obtained a limit-cycle solution, the analysis proceeds along very similar lines to that in the 
previous section. The perturbations in p and 7 are simple to derive and, since Eq. (|52|l is linear, they give exact 
solutions to the DWRG equation: 



1 



V3{p,j,A) 



V,^°\p,l,A) 



62,1, mA ^ p 7 . 

n.ni—O 



(58) 



The main features of the power-counting are unaffected by the addition of the new scale 7. The LO perturbation 
is still marginal, and leads to a family of limit-cycles parametrised by Co,o- Although it is more difficult to show, 
the analysis of the infinite-scattering-length case suggests that the marginal perturbation Co.o and the scale p, must 
represent the same three-body physics in the EFT. 

The previous power-counting is supplemented by terms proportional to powers of 7. In many systems, such as those 
of interest in nuclear physics, the two-body scattering lengths and hence 7 are fixed. This means that these terms 
can not be determined separately but must be absorbed into the coefficients of the energy-dependent perturbations. 
However in atomic systems with tunable Feshbach resonances, it may be possible to vary 7 and hence to disentangle 
the energy- and 7-dependent perturbations. 
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The full set of perturbations also includes ones that depend on the relative momenta. The forms of these have 
been derived and their eigenvalues show that they occur at the same orders as the corresponding en ergy -dependent 
perturbations. The full expressions are very unwieldy and so the interested reader is referred to Ref. \il\ . 

Having obtained a solution to the RG equation we can substitute it back into the LS equation to obtain the 
corresponding DW T-matrix. The algebra again follows that for systems with infinite scattering length. In the 
present case, the truncated Green's function involves an integral around the branch cut corresponding to elastic 
scattering of a single particle on a bound pair below the three-body threshold. This 2-1-1 continuum cut runs along 
the imaginary q-axis from ij, as shown in Fig. [3 Using the form for the three-body force given in Eq. ()58|l we find 
that the DW T-matrix element for 2-1-1 elastic scattering can be written 






-g2{p,l,P*)+ E C2„,V"7"- (59) 



p,i7K'3V-P/'l "^ p,i7/ n,m=0 

Using the properties of 02 in Appendix B, one can show that T^{p) respects unitarity for p^ < 0, as it should. 
For p^ > 0, the coupling to three-body breakup channels, mean that T-matrix for 2 -I- 1 scattering does not repect 
unitarity on its own. The same effective potential could be used to calculate amplitudes for breakup and 3 — > 3 
scattering. However, in practice this is likely to require general momentum-dependent potentials to describe all 
possible interaction channels. 

The DW T-matrix can be related to the elastic-scattering phase shift, S, by 

. = -— \cot(6{p) - 62ip)) - i], (60) 

where 62 is the phase shift produced by pairwise forces only. This allows us to rewrite Eq. H59|l in the form of a DW 
effective-range expansion, 

00 
^V2{p,j,p.)[cot{S{p)^S2{p))-i]-g2{p,l,P*)^- E ^2„,„^/"7"- (61) 



V. DISCUSSION 

In this paper we have applied a DW version of the RG to three-body scattering. This is based on cutting off the full 
solutions in the presence of long-range interactions Il4!|. In three-body systems, these long-range forces are generated 
by point-like interactions between pairs of particles. By imposing the cut-off on the DW's we ensure that it affects 
only the contributions of the short-range three-body forces, and does not alter the two-body physics. Demanding that 
observables are independent of the cut-off then leads to an RG equation for the renormalised three-body interactions. 

As shown by Efimov 26] , a three-body system with contact interactions can be described by an ISP in the region 
where they are all close together. Hence the RG in presence of this rather singular potential controls the behaviour 
of the three-body forces and determines their power counting. 

For systems of three bosons or two neutrons and a proton with total spin i, the resulting ISP is attractive. This 
leads to quite different behaviour compared with previous examples studied with this type of RG [H Il4| . The 
wave functions show oscillatory behaviour at short distances and a self-adjoint extension is needed to fix their phase 
and hence give well-defined eigenfunctions |29l 133. ISJL l35l . These oscillatory functions control the RG fiow of the 
short-distance potential, which tends to a limit cycle I38I rather than a nontrivial fixed point. 

These results are general for an attractive ISP in any number of dimensions. The lack of a ground state for such a 
potential means that it has deeply bound states that lie outside the domain of validity of our effective theory. This 
means that the cut-off should truncate not only high-energy scattering states but also these deeply bound states. This 
leads to a three-body force that has periodic discontinuities, but is single- valued. 

The resulting power counting for perturbations around the limit cycle starts with a marginal term. From its effect 
on observables, we have seen that this term is equivalent to changing the parameter in the self-adjoint extension or, in 
other words, the phase of the short-distance wave functions. To express the resulting power counting in a form which 
matches that usually used '7', 'o] , we can assign an order d = ^ — 1 to a term in the rescaled potential proportional 
to A*^. Then a term proportional to p^" (or E") appears at order d = 2n — 1. The coefficents in our potential have 
direct, and relatively simple, connection to scattering observables. 

The same behaviour is found in an attractive three-body system with finite two-body scattering length, 02, since the 
wave functions tend to the Efimov form at short distances (R <C 02). Hence the RG flow has a very similar pattern. 
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tending to a limit cycle. In this case the power counting should be extended to include terms involving powers of 
the low-energy scale 7 = l/a2- A term proportional to ^^"7™ is of order d — 2n + m — 1. For the energy-dependent 
perturbations, this counting agrees with that found by Bedaque et al. J23 from the STM equation. Because of the 
clean separation of the short- and long-range physics in the DWRG approach we are able to state the counting in an 
algebraically much simpler way. 

We are also able to extend the counting to terms involving powers of 7. This may not be relevant to applications 
in nuclear physics, where the two-body scattering lengths are fixed, but in atomic systems it may be possible to use 
Feshbach resonances to vary the scattering length, and hence to disentangle the two types of perturbation. 

The leading, marginal three-body force, or equivalently the choice of self-adjoint extension, is of interest since it 
shows that one piece of three-body physics is needed to determine low-energy three-body observables. As noted by 
Bedaque et al. [13, this provides a way to understand the old observation of the "Phillips line," a correlation between 
the J = ^ nd scattering lengths and triton binding energies predicted by model nucleon-nucleon potentials 5] . 

At leading order {d = —1), three-body observables are determined by the two-body scattering length and the 
marginal three-body force. The first correction to this arises from the two-body effective range, which appears at 
order d — 0. Including this has been shown to give good agreement with with J = -^ nd observables, in calculations 
using the STM equation 22] and an equivalent equation in coordinate-space |l7l |. The first energy-dependent three- 
body force appears only at next-to-next-to-leading order, d = 1 in this counting. 

It should be possible to extend our approach to three-body scattering above the breakup threshold. However this will 
introduce couplings between various channels and so will involve momentum-dependent perturbations in the effective 
potential. This is likely to require similar analyses to those which have been applied to simpler coupled-channel 
systems, such as that considered in Ref. |45j . 
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APPENDIX A: EFIMOV WAVE FUNCTIONS 

We summarise here the essential elements of Efimov's approach to the three-body problem 2d|, since this provides 
the wave functions we need to elucidate the scaling of the the three-body interactions. 

We consider the case of s-wave scattering. In general, the wavefunctions l^p.fc) niust be found using the Faddeev 
equations or equivalent (see, for example, Ref. [43j). In the Faddeev formalism, the wavefunction is broken into 
three components according to the pair of particles that interacted last. The component in which particles 2 and 
3 interacted last is denoted V'i ^(''23, ''i), where r23 and ri are the usual Jacobian coordinates for three bodies with 
equal masses. 

For particles interacting via pairwise contact interactions, Efimov observed that V'l*l.(''jfci^i) satisfies a free two- 
dimensional Schrodinger equation, subject to a boundary condition relating the wave function at the points Vjk — 
and ri — rjk/2 where pairs of particles interact. This boundary condition takes the form 



/,(») 



dtpi 



+ ^^«(i?,V3) = £^S(i?,0), (Al) 



where A is a factor related to the wave-function symmetry, and R and Lpi are the hyperradius and hyperangle respec- 
tively, and are defined by 



4 „ „ ^ VSr 



■jk 



i? = ^ -r2 + r% , (p, = arctan ^^ . (A2) 

For three bosons, the full wavefunction can be expressed as 

*M(^,^)-E]^2-^^i:i(^>^o, (A3) 



i=l 
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where H. represents five general hyperspherical coordinates that complement R. (The other coordinates could be ipi 
and two angles each for r^ and Tj^., but their exact specification will not be needed here.) 
We shall assume that the DW's are normalised by 

/•oo ^ 2 

/ R^dR dn^p^k{R,n)^p',k'{R,n) ^ —Sip''-p'^)d{k-k'), (A4) 

/•OO P 

/ R^dR dn'^p^,^{R,n)y5jp,^,^{R,n) ^ -s{p^-p'^), (as) 

/ R^dR I dn-^n{R,^)'i'n{RM) = 1- (A6) 

Because the boundary condition couples R and ipi, the equations are, in general, extremely difficult to solve. 
However when i? <C 02 the boundary condition separates. In this limit we may label the states in terms of the 
centre-of-mass energy and the hyperangular momentum, s. We can write the solutions in the separable form 

Tp^j^liR, ip^) = As sin (y - s(/7,) Us{p, R). (A7) 

This satisfies the angular Schrodinger equation subject to the boundary condition of vanishing amplitude at ipi = 7r/2 
(ri = 0). The boundary condition, Eq. HA1|I . now results in a transcendental equation for s, 

STT 8A . STT , . ^, 

scos — = — ;=sm — . (A8) 

This equation was also obtained much earlier by Danilov ,44] using the momentum-space equation derived by Skor- 
niakov and Ter-Martirosian |2^ . 

The resulting radial Schrodinger equation has the form 

d'^Us{p,R) 1 dus{p,R) s'^ I T,\ 2 [ T,\ /Am 
dR^ - R-^iR- + ]^"^(^'^) =^ -^(^'^)' (A9) 

and so contains an ISP whose strength is determined by s, and hence depends on the symmetry parameter. For 
three bosons or a neutron and deuteron with spin-1/2, the symmetry parameter is A = 1, while for three nucleons 
with spin-3/2 it is A = —1/2. In the latter case, Eq. I|A8|I has real solutions and inverse-square potential provides a 
repulsive "centrifugal barrier" . The power counting for short-range interactions in the presence of this potential can 
be derived using the RG method of Ref. 14| . 

The case with A = 1 is more interesting, since the lowest solutions to Eq. IJA8|I are imaginary. 



s = ±iso ==±a.006--- , (AlO) 

and so they correspond to an attractive ISP. For small R the solutions to the corresponding radial equation have the 
forms 

Us,{p,R)o,R^''\ (All) 

These are both equally (ir)regular as i? ^ 0. They can also describe flux disappearing or being created at the origin, 
corresponding to the classical "fall into the centre" which is possible for this potential JS^I- To obtain well-defined 
wave functions which respect ffirx conservation, we need to impose a boundary condition on the wave functions for 
i? — > 0, requiring equal admixtures of the two complex solutions and fixing their relative phase. Mathematically, this 
is known as choosing a self-adjoint extension of the Hamiltonian in Eq. HA9|I . as discussed by Bawin and Coon |29| 
and references therein. 

The full solutions to the Schrodinger equation with an attractive ISP can be written in terms of Bessel functions 
of imaginary order, 

uso {P,R)^ \l^^pr-TTT-r- 7W\ f ^'"^"^ -^-0 (pR) - ^-'"'(p) J_,,„ (pi?) 1 , (A12) 

V 2 2i\sin(r][p) + tnSo/2)\ I i 



where p = y/ME and the normalisation has been fixed by requiring that Usg {p, R) ~ sm{pR + 5)/\R for large R. For 
small R, these wave functions have the form 



Us. iP, R) - 1 / ^-^/^^^?^^^H:7TW ''<'" 1^^^^) + ''(P) - ^) ' (^13) 

y So[cosh(7rso) - cos[2r](p))\ 
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where 

6l = argr(l + iso) + Soln2. (A14) 

Demanding that these waves tend to a common, energy-independent form as i? ^ imphes that 77 (p) must be of the 
form 

'7(p) = -Soln(p/p,), (A15) 

where p, is a parameter which fixes the phase of the sine-log oscillations (or equivalently specifies the self-adjoint 
extension) . 

The physical meaning of p^ can be seen by noting that the S'-matrix for this system is 

^2u.i,)^^M'nip)-^-^on) (AI6) 

sm(77(p) -f inSo/2) 

This has a pole at p = ip^,, implying that E — ~p1/M corresponds to a bound state. In fact, as shown by Efimov, 
there is an infinite tower of bound states with energies E = —p\lM where 

P„-P*e"^/^«. (A17) 

The bound states accumulate at zero energy and extend downwards in a geometric pattern with no ground state. 
The wave functions of the bound states are 



ui:Hi?) - \r-^^^^p..K,s,xvnR\ (AI8) 

y ""So 
where Km{x) denotes a modified Bessel function of the third kind. Near the origin these have the form 



V2 

I H ] r^ 

So 



wloH^) ~ —Pn{-ir+^ sm{solnp,R~ 9) . (A19) 



s, 



APPENDIX B: RESCALED PROJECTION OPERATOR 

The short-distance Green's function defined in Eq. (|4()|l has the spectral representation 

72 p'' -cf + te [ 71" 7o 9 J 2 ^ p^ + pi 

Using the result, 

«^oJo xo-x±ie Jq Xf)-x 

where V denotes the principal value, we find that, for real p, 

i 2 f^ dk 

-[G2ip,l,P*)-G2i-p,l,P*)]^'D2{p,-/,p*) + - — X'3(p, fc,7,p*). (B3) 

n TT Jo P 

In this equation ^2(— P, 7,?*) is found by analytically continuing p through the upper half of the complex p plane 
to —p. Since, by its definition, G2ip,J,P*) is real for pure imaginary p, its value at negative real p is the complex 
conjugate of that at positive real p and hence G2i—p, 1,P*) corresponds to a — ie prescription at the propagator pole. 
By rescaling Eq. (|B3|I . we obtain an expression for the rescaled projection operator defined in Eq. (|5!^ : 

P(7,A) = ^(^2(7,A)-^2(-7,-A)), (B4) 

where Q2 is the rescaled Green's function defined in Eq. H49|l . 
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That rescaled Green's function has the representation 



G2{i,A) 



r die)- 4-- 'D2(l,qA)+^{f)- f dkvJ^l,qA 

J-r- l-q^+ie \q J tt Jq \q q 



'-7' 



The analytic properties of this function are needed since it forms part of the integrand H{q) of Eq. (|55|) . which is 
used to construct the RG hmit-cycle solution. From its spectral representation, we see that ^2(7, A) has poles at the 
bound state momenta, q = ipn{A). The integral term results in a branch cut down the imaginary axis from q = i-y 
and then along the positive real axis. 
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